Age and Dust Degeneracy for Starburst Galaxies Solved? 
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O : ABSTRACT 

o\ 

0\ • A spectral evolution model of galaxies that includes both stellar and dust effects 

^ . is newly built. xApplying the model to 22 nearby starburst galaxies, we have shown 

5^1 that far infrared luminosity of galaxies helps to break the age-dustiness degeneracy. We 

P \ have derived a unique solution of age and the dustiness for each starburst galaxy. The 

resulting starburst ages and optical depths are in the range 10 < f (Myr) < 500 and 
^ . 0.5 < Ty < 5.0, respectively. The result is robust and is almost independent of model 

assumptions such as dust distributions, extinction curves, and burst strengths. With the 
kp \ rapidly growing sensitivity of submillimeter detectors, it should become possible in the 

'j_j ■ near future to determine the age and Ty of star-forming galaxies at redshifts 2; ~ 3 and 

^ beyond. Accurate estimates of Ty for Lyman-break galaxies and high-z galaxies might 

require a substantial revision of the previously claimed picture of star formation history 
over the Hubble time. 

Subject headings: dust, extinction — galaxies: ISM — galaxies: starburst — galaxies: 
stellar content - radiative transfer 



1. INTRODUCTION 

Starburst (SB) galaxies in the nearby Universe are 
the best sites to study dust. Star formation in the SB 
galaxies should occur in localized regions, where mas- 
sive OB stars are deeply embedded in gas clouds that 
contain significant amounts of dust on average. The 
bulk of the photons emitted from burst populations 
are absorbed and partly scattered by dust, and are 
re-radiated in at far infrared (FIR) wavelengths. The 
intrinsic spectral energy distribution (SED) of the SB 
population is thus considerably modified in the UV. 
This applies not only to the nearby SB galaxies, but 
also to star-forming galaxies at high redshifts, since a 
timescale of heavy element synthesis in massive galax- 
ies at an initial burst of star formation should be as 
short as ~ 10^ yr (e.g., Arimoto & Yoshii 1987), and 
thus galaxies at high redshifts should contain as much 
dust as the nearby galaxies do. In recent years, there 
have been many attempts to derive the history of the 
star formation rate (SFR) in the Universe, e.g., by 
converting the UV light densities to the SFRs of the 
CFRS (Canada France Redshift Survey) galaxies at 
< z < 1 (Lilly et al. 1996) and the HDF (Hub- 
ble Deep Field) galaxies at 1 < z < 4 (Madau et al. 
1996; Connolly et al. 1997). However, without pre- 
cise estimates of the dust extinction and stellar age, 
it could be misleading to use the observed UV flux 
density as an indicator of the SFR, because the UV 
light is most sensitive to dust obscuration (Meurer 
et al. 1997; Rowan-Robinson ct al. 1997; Calzetti 
1998). Unfortunately, the UV-ncar infrared (NIR) 
SEDs are degenerate in stellar age and dust attenu- 
ation (Gordon, Calzetti, & Witt 1997). Complicated 
analyses of the SEDs are required to solve this degen- 
eracy by incorporating the dust effects into a model 
of galaxy spectral evolution. The nearby SB galax- 
ies are most suitable to calibrate the dust effects in 
actively star forming environments. Once the degen- 
eracy is solved, an attempt can be made immediately 
to determine the ages of Lyman-break galaxies (Stei- 
del et al. 1996) at high redshifts which are obviously 
too red to be dust-free star-forming galaxies (Meurer 
et al. 1997). It would also be possible to derive the 
age of primeval galaxy candidates in which the bulk of 
stars of present-day ellipticals formed in a relatively 
shortlived starburst for the first time (Hughes, Dun- 
lop, & Rawhngs 1997). 

Gordon et al. (1997) studied the UV-NIR SEDs 
of 30 SB galaxies which manifest star-forming activ- 



ity by the strong UV SEDs and intense line emis- 
sions. These arc irregular and spiral galaxies with 
disturbed morphologies locating at a median distance 
60 Mpc (hereafter Ho=50 km s~^ Mpc^^). Gordon 
et al. (1997) found that these galaxies show little 
evidence for a bump at 0.22/i?n, suggesting that the 
extinction curve similar to that of the Small Magel- 
lanic Cloud (SMC) is more suitable than that found 
for the Milky Way (MW). The best geometry that 
describes the SB region was found to have the inner 
dust-free cavity having star clusters surrounded by 
the dust shell with local clumpy dust condensations. 
Gordon et al. (1997) tried to get the ages for two SBs 
by using the UV-NIR colors, but failed to determine 
them uniquely. This suggests that the SB ages cannot 
be determined from the UV-NIR SEDs alone. 

The aim of this work is to investigate the SB galax- 
ies studied by Gordon et al. (1997) with additional 
IRAS data and to derive the accurate age and the dust 
attenuation of the SB population. In §2, we briefly de- 
scribe our model, in §3 wc describe our data sample, 
and in §4 we demonstrate that we are successful in 
disentangling the burst age and the dust attenuation 
in the UV-FIR SEDs, and derive the burst parame- 
ters for 22 individual SB galaxies. In §5, we discuss 
cosmological implications of our results and conclude 
this article in §6. 

2. Model of starburst galaxy 

Even in the SB galaxies, the burst population does 
not always dominate the galaxy's light, to which the 
underlying population contributes signiflcantly, and 
sometimes exceeds the burst light if the SB is already 
aged. We, therefore, adopt a reahstic model of star- 
forming spirals by superposing the burst population 
on the underlying galaxy. To model the SEDs of SB 
galaxies, both a stellar evolutionary synthesis model 
and a dust radiative transfer (RT) model arc needed. 
We use the stellar evolutionary synthesis model of 
Kodama & Arimoto (1997), the latest version of Ari- 
moto & Yoshii (1986; 1987). We calculate accurately 
the radiative transfer of the stellar photons through 
the dust by solving the transfer equations (Takagi, 
Vansevicius, & Arimoto 1998). Main ideas concern- 
ing the structure of SB galaxies are borrowed from 
Witt, Thronson & Capuano (1992) and Gordon et al. 
(1997). 



2.1. Structure of the starburst galaxy 

We have adopted a two-component model for the 
SB galaxy: the galaxy consists of the SB region and 
the underlying population and is characterized with 
a free parameter fsB , a- mass ratio of the SB to the 
total galaxy. The radiative transfer and the evolu- 
tionary models are calculated for both components 
independently. We assume that the densities of star 
and dust in the underlying population follow the mod- 
ified King's distribution: 
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where tq means a core radius of the galaxy. The 
distribution of stars in the SB region is modeled by a 
similar but steeper function: 
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where tsb means the core radius of the SB region. 
For the underlying and SB regions, the ratio of the 
total radius, Rtoti to core radius, Vcorei is fixed to be 
log{Rtot/rcore) — 2.2. The distribution of dust in the 
SB is assumed to be homogenous for the entire region. 
Both stellar and dust distributions are integrated up 
to the same limiting radius. 

As a test of validity of the geometry of our models, 
we have calculated a few models following the sugges- 
tion by Gordon et al. (1997). They assumed the SB 
model as the dusty sphere with the inner dust-free 
cavity where the stellar population is residing. No 
significant differences are found for relatively small 
optical depths tv <; 5, where ry is integrated from 
the center to Rtot- 

The model calculations are carried out for fss — 
0.1 — 1.0. In this paper, the main attention is drawn to 
the two cases Jsb — 0.23 and 1.0, which means 23% 
and 100% of galaxy mass are involved in the starburst. 
The cases with ]sb < 0.2 fail to cover the entire re- 
gion of the observed colors of the sampled galaxies in 
terms of age and ry. For the underlying, we adopt 
typical parameters for an Sc galaxy, i.e., a mass of 
1O"'^^M0, and a radius of 10 kpc. If we assume 30% 
SB mass to the underlying galaxy, it gives jsB = 0.23, 
the SB radius 6.6 kpc, and the mass 3 x lO^^M©, re- 
spectively, if the the same average density is adopted 



in the SB region as in the underlying galaxy. The to- 
tal mass of the galaxy is then 1.3 IO^^Mq. Such a size 
of the SB region is in agreement with the area cov- 
ered by the diaphragm for the observations presented 
in Gordon et al. (1997). 

The numerous tests with different geometrical pa- 
rameters reveal only small differences in the modeling 
results. However, in cases with ry > 5, the precise 
simultaneous determination of the age and dustiness 
of SB is model dependent, and have to be carefully 
calibrated. 

2.2. Radiative transfer in dusty galaxy 

The radiative transfer in the galaxy and SB is cal- 
culated using the RT code developed by Takagi et al. 
(1998). The code is based on an assumption of spher- 
ical symmetry with arbitrary radial distributions of 
stars and dust. Direct integration along the rays is 
adopted. The steps of the mesh (100 shells) are de- 
termined automatically, depending on the stellar and 
dust distributions, to fulfill the energy conservation. 
The ray tracing is done in the way suggested by Band 
& Grindlay (1985). The thermal emission from dust 
is calculated by assuming an equilibrium of dust par- 
ticles in the radiation field and by taking into account 
self-absorption. Up to six isotropic scattering terms 
are calculated, and losses of energy due to the scat- 
tering do not exceed 0.3%. The total energy losses in 
the present models are less than 1%. The isotropic 
scattering is assumed for simplicity without introduc- 
ing significant errors, since we analyze optically thick 
cases. We should note that the thermal emission, 
which is crucial for the present study, is independent 
of the adopted phase function. 

The results were tested by comparing with those 
obtained by Witt, Thronson, & Capuano (1992) and 
the output of RT benchmark by Ivezic et al. (1997), 
and were found to be consistent. Also no significant 
discrepancies arc found between our results and those 
obtained by DUSTY (Ivezic, Nenkova & Elitzur 1997) 
for the model of a star enshrouded in a dusty envelope 
with a central dust-free cavity. 

For the present investigation, Ty is adopted as a 
free parameter, i.e., we mimic the amount of dust by 
changing the optical depth of our models. That means 
the mass of dust depends on the details of model pa- 
rameters. This is why we do not derive directly the 
mass of dust but refer only Ty ■ Any extinction law 
can be applied in the code, but for the present inves- 



tigation we limit ourselves with the MW and SMC 
extinction curves (Pei 1992). The albedo is taken 
from Pei (1992) for both of the extinction curves. We 
truncate the albedo curve for A > 7/iTO, without in- 
troducing any significant error, in the calculation of 
the scattered light. 

2.3. Stellar evolutionary synthesis 

The evolutionary synthesis code developed by Ko- 
dama & Arimoto (1997) is used for the modehng of 
the luminosity and chemical evolution of the galaxy 
and SB. The code is based on the up to date ingre- 
dients, such as stellar isochrones, spectral libraries, 
and chemical yields of various elements, and has been 
successfully used in a numerous publications (e.g., 
Kodama et al. 1998; Kodama & Arimoto 1998). 
The code gives the evolution of synthesized spectra 
of galaxies (the wavelength range 0.09 — 3000/zm) in 
a consistent manner with the galaxy chemical evolu- 
tion, i.e., the effects of stellar metallicity are explicitly 
considered (Arimoto & Yoshii 1986). 

The underlying population of an Sc galaxy is mod- 
eled by following the prescription given by Arimoto, 
Yoshii & Takahara (1992), i.e., taking the values of all 
free parameters for an infall model; a SFR k — 0.190 
Gyr~^, an infall rate a = 0.148 Gyr^^, Salpeter ini- 
tial mass function (IMF) with a slope x ~ 1.35, and 
lower and upper stellar masses, m£ ~ O.IOM0 and 
niu = 6OM0, respectively. We adopt 15 Gyr as the 
age of galaxy because the observed sample consists of 
nearby galaxies. 

The SB population is modeled by assuming a 
closed box model with the SFR of ksB — 10.0 Gyr~^. 
The SFR adopted is equivalent to that of a giant el- 
liptical galaxy (Arimoto & Yoshii 1987). The neb- 
ular continuum is neglected in the present models 
since we limit ourselves to ages t > 10 Myr when 
the amount of ionizing photons (A < 0.09/i77i) is neg- 
ligible except for the very beginning of the bursts 
at i < 10 Myr (Leitherer & Heckman 1995; Fioc 
& Rocca-Volmerange 1997). The IMF for the SB 
is assumed to be Salpeter-like with x = 1.10, and 
me = O.IOM© and m„ = 6OM0, respectively. Such 
a flat IMF has been suggested for spheroidal galax- 
ies (e.g., Arimoto & Yoshii 1987; Kodama & Arimoto 
1997) and for the OB associations in the Milky Way 
(Massey, Johnson, & DeGioia-Eastwood 1995). 

The metallicities given by Calzetti, Kinney, & 
Storchi-Bergmann (1994) for the SB galaxies are at- 



tributed to the ionized gas, and the averaged value for 
the sampled galaxies {\og{0/H) -H 12 = 8.7 ± 0.5) is 
very close to solar. Nevertheless, we calculate our 
models with the initial metallicity, Zi = 0.02 and 
0.008, there is no significant differences. Therefore, 
the initial metallicity of the SB region is adopted to 
be solar. 

The SEDs of the underlying galaxy and the SB 
region at different ages are shown in Fig.l. The thick 
solid line represents the underlying at 15 Gyr. For the 
SBs, all SEDs are plotted for Ty — 1.0. The optical 
depth for the underlying galaxy is fixed at Ty = 0.5 
throughout the paper. 

3. Observational data 

The observational data of SB galaxies are taken 
from Gordon et al. (1997). They presented the ho- 
mogeneous sample of aperture-matched observations 
of 30 SB galaxies from the far-UV to NIR. The size of 
the diaphragm of the observations is 10 x 20 arcsec, 
which is equal to ~ 4.5 kpc at their sample's median 
distance of 60 Mpc. 

We have supplemented the NIR colors for 4 galax- 
ies (NGC1510, NGC5236, NGC5253 and NGC7552) 
by taking data from available multiaperture photom- 
etry (de Vaucouleurs & Longo 1988), which allows 
us to reduce the observed magnitudes to the effective 
aperture used by Gordon et al.(1997). 

The IRAS flux densities are taken from NED (NASA/IPAC 
Extragalactic Database) and have been color cor- 
rected. Two galaxies (Mrk 309 and Mrk 542) in the 
sample have only upper limits for the lOOfim flux den- 
sity. It is important to note that the IRAS colors are 
derived for the total galaxy while our sample (UV- 
NIR) is aperture limited. Nevertheless, taking into 
account the conclusion drawn by Calzetti et al.(1995) 
concerning the strong central concentration of FIR in 
the galaxies of the present sample, we believe that the 
necessary aperture correction does not change signif- 
icantly our main conclusions. 

Finally we end up with a sample complete in B, 
H, and IRAS 60 and 100 /im fluxes for 22 galaxies. 
The choice of present colors B, H, and FIR is done 
after an analysis of all the observational data given 
by Gordon et al.(1997). Only in the case of B-H and 
Lfir/ Ln, the disentangling of the age and dustiness 
is successful. LpiR is the FIR luminosity derived from 
60 and 100 /iw fluxes (Londsdale et al. 1985). Table 
1 gives the observations and derived parameters of 



the SB galaxies: column (1) gives the name of galaxy, 
columns (2) and (3) indicate B-H and log(LFiR,/iH), 
respectively, and columns (4)-(ll) give the derived 
age in Myr and the optical depth (see §4) for both 
the MW and SMC extinction curves. 

4. Results 

The comparison of model results with the observed 
properties of SB galaxies is divided into two steps: a 
determination of the age and dustiness of SB galaxies, 
and a comparison of the observed SEDs with model 
predictions for the age and the opacity sequences. 

4.1. Simultaneous determination of the star- 
burst age and dustiness 

The calibration of SB galaxies with fsB — 0.23 by 
the age and dustiness for the MW extinction curve 
is shown in Fig. 2a. The data points are taken from 
Table 1. The large circle indicates the underlying 
galaxy. An area indicated by a dashed curve repre- 
sents the region occupied by Coma spirals (Gavazzi, 
Randone & Branchini 1995). The solid lines show the 
isochrones with increasing ry from the left (ry — 0.5) 
to the right (ry = 10). The dot-dashed lines show the 
isoopaques with increasing age of the SB from the top 
(t = 10 Myr) to the bottom (i = 1 Gyr). 

It is obvious that FIR luminosity of galaxies helps 
to break the age-dustiness degeneracy, and we have the 
unique solution for the SBs with Ty < 5. The range 
of Ty validity depends on fsB assumed. Fig. 2b shows 
the same diagram as Fig. 2a but for fss = 1, i-e., 
the galaxy has 100% SB population. As we can see 
in Fig. 2b, the applicable range of ry is significantly 
larger, but the shape of the grid remains very similar 
to the case shown in Fig. 2a. Therefore, we expect that 
the assumed value of the free parameter fsB = 0.23 
does not infiuence our conclusions significantly. The 
derived ages and optical depths Ty for individual SB 
galaxies are given in Table 1 for fss = 0.23 and 1.00. 
We note that the burst age and the opacity are nearly 
the same for different values of fss- The resulting 
age ranges from t = 10 Myr to i = 500 Myr, which 
agrees well with Gordon et al.'s (1997) estimate for 
the present sample of SBs. 

Figs. 2c, d are the same as Figs. 2a, b, respectively, 
but for the SMC extinction curve. We note that the 
burst age and the opacity are nearly the same for dif- 
ferent extinction laws. We therefore can conclude that 
the present method is robust and does not depend on 



the details of extinction curve and the burst strength. 

4.2. SEDs of dusty starburst galsLxies 

The parameters given in Table 1 are used for group- 
ing of the SB galaxies by the age and the opacity. 
Two groups of galaxies are defined: (a) those with 
approximately the same ages {t ~ 20 — 30 Myr) but 
different optical depths {tv = 0.5 — 5.0), and (b) those 
with similar optical depths {ry = 0.9 — 1.5) but dif- 
ferent ages (i = 10 - 400 Myr). The SEDs inside 
these groups are averaged and plotted in Figs.3a-h. 
The model SEDs are superposed on the observational 
data without any fitting, but both are normalized at 
the H band to 1.0. 

Figs.3a-d and Figs.3e-h represent the two different 
extinction curves (MW and SMC, respectively). The 
SEDs in Figs.3a,b,e,f represent Ty sequences and in 
Figs.3c,d,g,h age sequences. The infiuence of the ex- 
tinction law is significant in both the UV and FIR 
spectral regions. The behavior of the SED in UV can 
be predicted from the shape of the extinction curve 
only in the case of fss = 1.0. A strong nonlinearity is 
observed for Ty > 2 in case of fsB — 0.23, where the 
stronger infiuence of fsB than of the extinction law on 
the UV part of SED is observed, since the significant 
obscuration in the UV hides SB population and the 
underlying galaxy, which is characterized by a MW 
extinction curve, becomes dominant. FIR emission is 
affected more regularly by the extinction curve and 
the burst strength, i.e., the FIR part of SED becomes 
cooler for the SMC extinction law and brighter for 
fsB = 1- As we can judge from the far-UV end of 
the SEDs, the extinction law is steeper than MW but 
flatter than SMC. Therefore, Calzetti et al. (1994) 
and Gordon et al. (1997) conclusions concerning the 
far-UV shape of the extinction curve in SB galaxies 
are confirmed. 

It is important to note that the observed difference 
in the UV part of SED is much smaller for the age 
sequence than for the Ty sequence, but the trend is 
opposite in the FIR region. Along the Ty sequence, 
the UV-NIR part of SED becomes fainter, while the 
FIR part becomes more luminous. On the other hand, 
both the UV-NIR and FIR parts of SED becomes 
fainter along the age sequence as the stellar popula- 
tion gets older. 

Figs.3a-h show that B,II and FIR colors in gen- 
eral give a reasonable estimation for the SEDs of SB 
galaxies from the UV to FIR. Nevertheless, we have 



to stress that Figs.3a-h are given without any fitting 
of the parameters, i.e., they represent a general be- 
havior of our models used to determine the age and 
TV- Moreover, the observational data are averaged 
for different age and dustiness groups. The fine tun- 
ing of model parameters for each particular galaxy is 
a subject for a forthcoming paper. 

5. Discussion 

Gordon et al. (1997) suggested that the dust dis- 
tribution in SB is best described by a clumpy me- 
dia, which is not adopted in our modeling. However, 
we should note that our Figs.2a-d would not change 
even if we take into account clumpy media, because 
they can be scaled in B-H and Lfir/Lh indepen- 
dently of the dust attenuation; thus our arguments 
are safe so far as we do not discuss the mass of dust 
which is model dependent. Moreover, Witt & Gordon 
(1996) suggested that the high conversion efficiency of 
UV/optical radiation into far-IR thermal dust emis- 
sion observed in SB galaxies and mergers is due in 
large part to a change in the structure of the inter- 
stellar medium in such systems toward a more ho- 
mogeneous widely extended distribution. Taking into 
account this suggestion and the curves correspond- 
ing to homogeneously extended distribution in their 
Fig. 19 (the density ratio of interclump-to-clump me- 
dia, k2/kl > 0.3) , we can expect that our definition 
of TV for homogeneous distribution of dust is linearly 
related with their tv calculated for the clumpy media. 

Figs.3a-h clearly reveal two problems: a discrep- 
ancy of the model SEDs and the observed ones in UV 
region, and a difference of the predicted and observed 
fluxes in mid-infrared region (MIR). 

The first point was discussed extensively by Gor- 
don et al. (1997). We have confirmed their conclusion 
that the extinction curve in the SB galaxies is steeper 
than in MW. Therefore, from the comparison of the 
ages and ry's in Table 1, we conclude that the uncer- 
tainty of the extinction curve in the UV is not serious 
for the present study. However, the deviation of the 
MW and SMC extinction curves adopted by us from 
the real one cannot significantly change our results, 
because the main diagrams (Figs.2a-d) are based on 
the bluest band B which is not influenced yet by the 
UV steepening of the extinction curve, and because 
the supplementary amount of the energy absorbed in 
far-UV does not increase significantly the FIR flux 
(Trewhella, Davies & Disney 1997). 



The second point, probably, can be explained by 
the shortcomings of our RT model, because the cir- 
cumstcllar dust and emission from very small grains 
are not included yet. The circumstellar dust mainly 
influences 12/iTO flux and is dependent on the num- 
ber of stars with dusty shells in galaxy. The influence 
of small grains on the IRAS colors was discussed by 
Boulanger et al. (1988). They concluded that IRAS 
12 and 25/Km colors depend on the non-equilibrium 
emission from small grains while 60 and lOO/im colors 
depend on the equilibrium emission from large grains. 
Extensive discussion on the applicability of so called 
"classical" grain model to flt the MIR (5 — 40/iTn) 
fluxes from the interstellar medium was given by 
Hclou, Ryter & Soifer (1991). Therefore, our model 
represents correctly the FIR emission (60 and lOO/xrn) 
from the grains in thermal equilibrium, which is cal- 
culated explicitly by our RT code. Moreover, we are 
using only the FIR colors for age-dustiness determi- 
nation, and can neglect the discrepancies between the 
predicted SEDs and the observed colors in the MIR 
region. 

Figs.3a-h indicate that the observed FIR flux ra- 
tios Sqq/ SioQ (dust temperature indicator) are smaller 
than the model predictions. This ratio depends on the 
dust temperature and on the emissivity law of dust 
grains. The dust temperature depends on the distri- 
bution of stellar density, which is rather approximate 
in our model, a flatter one would result in the cooler 
dust temperature. We assume the emissivity with the 
power /? = 2 (A"'', Rowan- Robinson, 1992) If /3 < 2 is 
taken, the discrepancy of the flux ratio would reduce 
(Hughes et al. 1997). However, such a fine tuning of 
the model is beyond a scope of this article, since both 
effects would not change significantly the total FIR 
flux used for breaking the age-dustiness degeneracy. 

We have demonstrated that the use of SEDs from 
far-UV to FIR is quite effective in disentangling the 
age and opacity of the SB galaxies. For deriving the 
opacity, Calzetti, Kinney, & Storchi-Bergmann (1996) 
and Calzetti (1997) assumed that the intrinsic line 
ratios arc given by a simple model with the optically 
thick limit, the so-called Case B of Ostcrbrock (1989), 
and estimated the opacity from the ratios of attenu- 
ated hydrogen emission lines. Although the opacities 
derived by two methods are correlated, the opacity 
given by the line method is smaller than the FIR 
method. This in turn should give systematically older 
ages for the line method. 

One should keep in mind that some of the SB 



galaxies assigned for ages younger than 20 Myr could 
well be much younger if their colors are significantly 
reddened by possible nebular continuum (Fioc & 
Rocca-Volmerange 1997). However, even if this is 
the case, the probability for galaxies younger than 10 
Myr to reside in this region should be very small, be- 
cause of short timescale for that stage. Additionally, 
if galaxies are UV-selected as in Gordon et al. (1997), 
it would be almost impossible to find such very young 
SB galaxies due to strong dust obscuration. 

Before we conclude this article, we briefly discuss 
two cosmological implications of the present results. 
First, if we define the primeval galaxies as massive 
galaxies at the onset of the first episode of major star 
formation, presumably progenitors of the present-day 
ellipticals at high redshifts, these galaxies should look 
like the SB galaxies in which burst population domi- 
nates the galactic light. Thus, it should be possible to 
determine the age and optical depth of the primeval 
galaxy candidates with a help of the rest-frame B- 
H vs. LpiYi/Ln diagram. These galaxies should be 
spectacular in the rest-frame FIR regime. At the 
highest redshifts, the FIR 60 and lOO/im fluxes conve- 
niently fall on into the region of submillimeter atmo- 
spheric windows at 350, 450, 750 and 850/i?7i. There- 
fore, a search for the primeval galaxies with infrared 
satellites and submillimeter telescopes, such as IRIS 
(Infrared Imaging Surveyor; Murakami et al. 1998) 
and LMSA (Large Millimeter Submillimeter Array; 
Ishiguro et al. 1994), would be promising. Indeed, 
the sensitivity of submillimeter arrays is expected to 
achieve the greatest advances in the immediate future 
(Hughes et al. 1997). 

Second, once the rest-frame B-H and Lpin/Ln 
are observed for high redshift galaxies including the 
Lyman- break galaxies (Steidel et al. 1996), radio 
galaxies (Dunlop et al. 1996), and QSO host galaxies, 
the ages and optical depths of these galaxies can be 
determined uniquely. Then, the resulting age should 
play the role of a cosmological clock that gives the age- 
redshift relation of galaxies at high redshifts and will 
constrain the cosmological picture of the Universe. 
The resulting ry should in turn give a correction fac- 
tor for the UV light, which will make it possible to 
derive the correct SFR at high redshifts. In other 
words, it could be misleading to derive the SFR for 
high redshift galaxies from the UV light alone. 



6. Conclusions 

In an attempt to break the age-dustiness degener- 
acy of the SEDs of starburst galaxies, we have built a 
new spectral evolution model of galaxy that included 
both stellar and dust effects. The evolutionary code 
of Kodama & Arimoto (1997) is used to model the 
spectral and chemical evolution of galaxies, and the 
radiative transfer is calculated using the RT code de- 
veloped by Takagi et al. (1998). Applying the model 
to 22 nearby SB galaxies sampled by Gordon et al. 
(1997), we have shown that the unique solution can 
be sought if the FIR luminosity is used for the star- 
bursts with Ty < 5. Except for two opaque galaxies 
with ry > 5, the resulting SB ages and opacities are 
in the range 10 < t (Myr) < 500 and 0.5 < tv < 5.0, 
when the burst mass fraction fsB — 0.23 is assumed. 
The result is robust and almost independent of the 
details of the dust distribution, the extinction curve, 
and the burst strength. Even in a case of the SB 
dominant galaxies, the age and Ty are uniquely as- 
signed on the rest-frame B-H vs. Lfir/^h diagram. 
For high redshift galaxies, these wavelength regions 
conveniently fall on into the NIR, MIR, and submil- 
limeter regions. With the rapidly growing sensitivity 
of the submillimeter detectors, it should soon become 
possible to determine the age and ry of star-forming 
galaxies at redshifts z ~ 3 and beyond. Accurate 
estimates of tv for the high redshift galaxies might 
require a substantial revision of the broad picture of 
star formation history over the Hubble time proposed 
by Madau et al. (1996). 
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Table 1 
Observed colors, derived ages and optical depths of SB galaxies 
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Note. — f The ratio of the SB mass to the total mass of galaxy, fsB 

: The ages and optical depths derived by the upper limit of LpiR. 
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Fig. 1. — The spectral energy distributions (SEDs) of SB region with a mass, 3 x lO^^M© , ry = 1.0, and different 
ages (thin Hnes) and of the underlying galaxy with a mass, IO^^Mq, tv = 0.5, and i = 15 Gyr (thick solid line). 
The ages of SB region are shown in the figure. The MW extinction curve is applied. 
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Fig. 2. — (a) The calibration of SB galaxies by the age and dustiness. All models are calculated with fsB — 0.23 
and the MW extinction curve. The filled circles are data points given in Table 1. The solid lines are isochrones 
with Tv = 0.5,1.0,2.0,3.0,5.0,7.0, and 10.0 from the left to the right. The dot-dashed lines are isoopaques with 
t = 10, 30, 100, 300Myr and 1.0 Gyr from the top to the bottom. 
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Fig. 2. — (b) The same as (a) but for fss = 1-0. 
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Fig. 2. — (c) The same as (a) but with the SMC extinction curve. 
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Fig. 2. — (d) The same as (c) but for fsB = 1-0. 
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Fig. 3. — Comparisons of the model SEDs with the averaged empirical SEDs within the groups of galaxies (see 
below for the definition of groups). All SEDs are normalized at H band, (a) The group of SB galaxies with similar 
ages {t ^ 20 r^ 30 Myr) but different opacities. The galaxies are divided into four sub-groups by the optical 
depths, TV — 0.5 (triangles), l.O(crosses), 3.0(circles), 5.0(diamonds). The theoretical loci are calculated with the 
age i = 30 Myr, fsB = 0.23, the MW extinction curve, and the optical depths ry = 0.5 (solid line), 1.0 (dotted 
line), 2.0 (dashed line), 3.0 (dot-dashed line), and 5.0 (long dashed line), respectively, (b) The same as (a) but for 
fsB = l-O- (c) The group of SB galaxies with similar optical depths (ry ~ 0.8 ^ 1.5) but different ages, i 2± 10 
Myr (triangles), 30 Myr (crosses), 100 Myr (circles), and 400 Myr (diamonds). The theoretical loci are calculated 
with TV — 1.0, fsB = 0.23, the MW extinction curve, and the ages 10 Myr (sohd line), 30 Myr (dotted line), 100 
Myr (dashed line), 300 Myr (dot-dashed line), and 1 Gyr (long dashed line), respectively, (d) The same as (c) but 
for fsB = 1-0. (e)-(h) The same as (a)-(d), respectively, but with the SMC extinction curve. 
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